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Abstract 



Three CP-violating phases of the 3x3 lepton flavor mixing matrix V are 
entangled with one another in the renormalization-group evolution from the 
seesaw scale (Ass ~ 10 14 GeV) to the electroweak scale (Aew ~ 10 2 GeV). 
Concerning the Dirac phase 5, we show that 5 = 90° at Aew can be radiatively 
generated from 5 = 0° at Ass i n the minimal super symmetric standard model, 
if three neutrino masses are nearly degenerate. As for the Majorana phases p 
and a, it is also possible to radiatively generate p = 90° or a = 90° at Aew 
from p = 0° or a = 0° at Ass- The one-loop renormalization-group equations 
for the Jarlskog invariant and two off-diagonal asymmetries of V are derived, 
and their running behaviors from Ass to Aew are numerically illustrated. 

PACS number(s): 14.60.Pq, 13.10.+q, 25.30.Pt 
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I. INTRODUCTION 



Recent solar [1], atmospheric [2], reactor (KamLAND [3] and CHOOZ [4]) and accelerator 
(K2K [5]) neutrino oscillation experiments have provided us with very robust evidence that 
neutrinos are massive and lepton flavors are mixed. The phenomenon of lepton flavor mixing 
can be described by a 3 x 3 unitary matrix V, commonly referred to as the Maki-Nakagawa- 
Sakata (MNS) matrix [6]. A useful parametrization of V reads [7]: 



V 



( c 12 c 13 s 12 c 13 S 13 \ 

— C 12 S 23 S 13 — S 12 C 23 e — S 12 S 23 S 13 + C 12 C 23 e S 23 C 13 

V — C 12 C 23 S 13 + S 12 S 23 e — S 12 C 23 S 13 — C 12 S 23 e C 23 C 13 / 



V o 




(i) 



where c^- = cos% and = sin% (for ij = 12,23 and 13). Note that the CP- violating 
phase 5 governs the strength of CP or T violation in normal neutrino oscillations, and it does 
not appear in the expression of (m) ee — the effective mass term of the neutrinoless double- 
beta decay 1 . On the other hand, the CP-violating phases p and er only affect (m) ee , and 
they can be rotated away if the massive neutrinos are Dirac particles. Current experimental 
data indicate 9 12 ~ 33°, 9 23 ~ 45° and 9 13 < 10° [8], but the phase parameters 5, p and 
o are entirely unrestricted. A variety of new neutrino experiments are underway, not only 
to measure the smallest mixing angle #13 and the Dirac phase 5, but also to constrain the 
Majorana phases p and er. 

While neutrino masses and lepton flavor mixing parameters can be measured at low- 
energy scales, their origin is most likely to depend on some unspecified interactions at a 
superhigh energy scale. For instance, the existence of very heavy right-handed neutrinos 
and lepton number violation may naturally explain the smallness of left-handed neutrino 
masses via the famous seesaw mechanism [9] at the scale Ass ~ 10 14 GeV. Below this seesaw 
scale, the effective Lagrangian for lepton Yukawa interactions can be written as 

-C = ElHWn - \e~ l H 2 ■ k ■ H?E C L + h.c. (2) 

in the minimal supersymmetric standard model (MSSM) 2 , where E L denotes the leptonic 
SU(2) L doublets, H 1 and H 2 are the Higgs fields, l R denotes the right-handed charged 
leptons, H2 = ia 2 H2 and El = ia 2 CE L with C being the Dirac charge-conjugate matrix. 
After spontaneous gauge symmetry breaking at the electroweak scale A E w ~ 10 2 GeV, 
we arrive at the charged lepton mass matrix Mi = vYi cos f3 and the effective Majorana 
neutrino mass matrix M v = v 2 Ksin 2 [3, where v ~ 174 GeV and tan/3 is the ratio of the 
vacuum expectation values of H 2 and Hi in the MSSM. The lepton flavor mixing matrix 
V arises from the mismatch between the diagonalization of Y\ (or Mi) and that of k (or 



1 Namely, (m) ee = \m 1 c 2 2 c 2 3 e 2tp + r m 2 s 2 2 c 2 3 e 2 ' la + m 3 s 2 3 \ is independent of 5. 

2 For the sake of simplicity, we assume the supersymmetry breaking scale Asusy to be close to the 
electroweak scale Aew- Even if Asusy/Aew ~ 1A the relevant RGE running effects between these 
two scales are negligibly small for the physics under consideration. 
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M v ). In the flavor basis where Y\ is real and diagonal, V directly links the neutrino mass 
eigenstates u 2 , f 3 ) to the neutrino flavor eigenstates (z/ e , z/^, v T ). The physical parameters 
at A S s an d A EW are related by the renormalization group equations (RGEs) [10]. It has been 
shown that the RGE evolution between these two scales may have significant effects on the 
mixing angle 6*12 and the CP- violating phases 5, p and a [11-15], in particular if the masses 
of three light neutrinos are nearly degenerate. It has also been noticed by Casas et al [12] 
and Antusch et al [13] that a CP-violating phase can be radiatively generated due to the 
RGE running from A S s to A EW . 

The reason for the radiative generation of a CP- violating phase is simply that three phases 
of V are entangled with one another in the RGEs. In other words, the running behavior 
of 5 depends on a non-linear function of 5, p and a. It is therefore possible to generate a 
non-zero value of S at Aew even if 5 = holds at Ass, provided the initial values of p and 
a are not vanishing. This observation opens a new and interesting window to understand 
possible connection between the phenomena of CP violation at low- and high-energy scales; 
e.g., the phase parameter governing the strength of CP violation in a long-baseline neutrino 
oscillation experiment could be radiatively generated from those CP-violating phases which 
control the leptogenesis of right-handed neutrinos [16] at the seesaw scale. 

The main purpose of this paper is to analyze the radiative generation of three CP- 
violating phases (<5, p, a) via the one-loop RGE running effects from A S s to A EW . Our work 
is different from the one done in Refs. [12] and [13] at least in the following aspects: 

• The phase convention of V taken by us in Eq. (1) forbids the phase parameter 5 to 
appear in the effective mass of the neutrinoless double-beta decay (m) ee . Hence it 
makes sense to refer to 5 as the Dirac phase. In contrast, the so-called "Dirac" phase 
in the "standard" parametrization of V used by Casas et al [12] and Antusch et al [13] 
enters the expression of (m) ee and is not purely of the Dirac nature. Our RGEs for 
three CP-violating phases turn out to be remarkably different from theirs. 

• We focus on the possibilities to generate S = 90°, p = 90° or a = 90° at Aew from 
5 = 0°, p = 0° or a = 0° at A ss . While 5 = 90° might imply "Maximal" CP violation 
in some sense [17], p = 90° or a = 90° will lead to a kind of large cancellation in 
(m) ee . The RGE running of 5 is of particular interest, because it means the radiative 
generation of leptonic unitarity triangles in the complex plane; i.e., three overlapped 
lines (sides) at Ass can evolve into a triangle at Aew- 

• The one-loop RGEs for the Jarlskog parameter (denoted by J) [18] and two off- 
diagonal asymmetries of V (denoted as A L and «4 R ) [19] are derived, and their run- 
ning behaviors from A S s to A E w are illustrated. These rephasing-invariant quantities 
measure the strength of Dirac-type CP violation and the geometrical structure of V, 
respectively. For example, A h = (or A R = 0) would imply V to be symmetric about 
its V^i-V^-Vrt (or V^-V^-Ki) axis. 

Because the radiative corrections to V are expected to be insignificant in the case that three 
light neutrinos have a strong mass hierarchy and tan (3 takes small or mild values, we shall 
mainly consider the quasi-degenerate neutrino mass spectrum in our numerical calculations. 

Section II is devoted to the one-loop RGEs for three neutrino masses, three lepton flavor 
mixing angles, three CP-violating phases and three rephasing-invariant quantities of V . A 
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detailed numerical analysis of radiative corrections to the Dirac and Majorana phases is 
presented in section III, where the RGE evolution of J , A^ and A^ is also analyzed. We 
give a brief summary of our results with some concluding remarks in Section IV. 



II. RGES FOR PHYSICAL PARAMETERS 



Below the seesaw scale, the effective neutrino coupling matrix k obeys the following 
one-loop RGE in the MSSM: 



lQn 2 ^ = aK+(Y l Y?)K + K(Y l Y?y 



dt 



with t = In (/x/ Ags) an d 



« = ~9 2 i ~ 692 + 6 (vl + y 2 c +Vt) , 



(3) 



(4) 



where g l and g 2 denote the gauge couplings, and (for f — u,c, t) stand for the Yukawa 
couplings of up-type quarks. In the flavor basis where Yi is real and diagonal, we have 
k = VkV t with k = Diag{/t 1? k 2 ,k 3 }. The neutrino masses at A EW read as m i = v 2 ^ sin 2 j3 
(for i = 1,2,3). It is then possible to derive the RGEs for (k v k 2 , /%), (#12, #23, #13) and 
(8, p, a) from Eq. (3), similar to the work done in Refs. [12-15]. 

For simplicity, we neglect the small contributions of y e and y^ to Eq. (3) in our calcula- 
tions. The RGEs of K i (for i — 1,2, 3) turn out to be 



6.K 1 K-y 

~dt ~ 16^ 



+ 2?/ T ^S 12 S 2 3 2c < 5C 12 C23S 12 S 2 3S 1 3 + C 12 C 23 S l:: 



dKo 



dt 16tt 2 



) 

a + 2y T (c 12 s 23 + 2c < 5C 12 c 23 s 12 s 23 s 13 + c 23 s 12 s 13 ^ 



dt 16tt 2 



a + 2 2/r C 23 c2 



t u 23 u 13 



(5) 



where c s = cos 5. The RGEs of % (for ij = 12, 23, 13) are found to be 



d6 



12 



y; 



dt 



16tT 2 { Cl2 

+Cl2 S (p-(T) 



( 2 _ 2 2 

C (p-cr) C 12 S 12 ^ S 23 C 23 S 13y 



C 23 S 23 S 13 



( 2 _ 2 2 

S (p-o-) C 12 S 12 \ S 23 C 23 S 13, 



L „ , f,2 _ „2 2 \ / _ 2 \ 

6 13 J \ C (8+ P -a) c l2 C (5-p+CT) tS 12j 

) ~~ ( S (5+p-cr) c 12 + S (5-p+cr) s 12) C 23 S 23 S 13 

( c (<5-p) s 12 s 23 ~~ C p C 12 c 23 S ls) — Cl3 S p ( s (5-p) s 12 s 23 + S p C 12 c 23 S 

~T~ ( C ((5-cr) C 12 S 23 + C cr C 23 S 12 S 13) — C23 S cr ( S (5-cr) C 12 S 23 ~ S cr C 23 S 12 S 13) 



C13 



C 23 S 12 S 13 



C 12 C 23 S 13 f ' 



d^23 = Vr^S 
dt 167T 2 



'(8-p) 



>13 



(< 



C (<5-p) S 12 s 23 C p C 12 c 23 S 13) + Cl3 S (<5-p) ( s (<5-p) s 12 s 23 + S p C 12 c 23 S 13y 



'12 



+ 



\5-a) 



s23 



( C (<5-ct) c 12 s 23 + C cr C 23 S 12 s 13) + C23 S (<5-<r) ( s (<5-ct) c 12 s 23 S ct C 23 S 12 s 13) C 12 j ) 
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d#i 3 2/rC; 



2:-! 



dt 16tt 2 

c, 



Cl3 



c (5-p) s 12 s 23 C p C 12 c 23 S 13 



) - Cl3*p ( 



s (<5-p) s 12 s 23 + S p C 12 c 23 S 13 , 



c 12 c 13 



+ 



£ ( C (5-cr) C 12 S 23 + C <t C 23 S 12 S 13) C23 S cr ( S (5-ct) C 12 S 23 S <r C 23 S 12 S 13) 



C 13 S 12f , (6) 



in which ^ - = (/^ — -) / + «■) (for = 12,23,13), c A = cos A and s A = sin A (for 
A = 5, p, a, 5 — p, 5 — a, p — a, 5 + p — a, 5 — p + a). In addition, the RGEs of 5, p and 
a are obtained as follows: 



dt levrM C12 



C (p-°") ( S 23 C 23 S 13) ( 



C (<5+p-<r) C 12 ~~ C (<5-p+cr) S 12 



2 \ C 23 S 23 S 13 



12 a 12 



Cl2 C (p-cr) 



f 2 2 2 \ f 2 , 2 

S (P-^) V s23 ~~ C 23 S 13j ~~ \ S (<5+p-<r) C 12 + S (6-p+a) S 12 



2 \ C 23 S 23 S 13 



r 1 C ( 5 -P) S 12 S 23 C p C 12 C 23 S 13 ) 
S13 

+C13 ( S (5-p) S 12 S 23 + S P C 12 C 23 S 13) 

~~ T~~ ( c (<5-o-)C 12 S 2 3 + C .C 2 3S 12 Si3 1 
1,23 x 7 



S p C 23 



C 12 S 12 
Si 



c 12 s 13 



~C23 { S (S-a 
d P _ Vr \ S (p-<?) 



) C 12 S 23 S o-C 23 S 12 S 13 



C p C 23 
c 12 s 13 
S o- C 23 
S 12 S 13 



\ C <r C 23 / 2 _ 



/ 2 2 2 \ _ 6 12 (2 2 \ 

\ C 12 S 12 S 13) S {S-p) „ l C 23 S 23i 

6 23 

/ 2 2 2 \ , S 12 / 2 2 \ 

\^ c 12 s 12 s 13j > c (<5-p)~ \^ c 23 S 23j 

/ 2 2 2 \ | C 12 (2 2 \ 

\ s 12 C 12 S 13,J + s (5-ct) l^ c 23 s 2Z) 



s 12 s 13 



2 2 \ 
C 12 S 13 j ~ C (S-a) 



u 12 
s 23 



( C 23 S 2s) 



( 



i 2 2 2 \ / 2 , 2 

^(p—cr) ^12^12 1^23 ^23^13 J I CfX-i-n—*-} Ci 9 Cf^^j.v* Si 9 ) Cq-jSq^S 







dt 16tt 2 [ C12 

— Cl2 C (p-cr) S (p-cr) C 12 S 12 ( S 23 — C 23 S 13) ~~ ( S (<5+p-cr) C 12 + S (8-p+o) S l2) C 23 S 23 S 



(<5+p-cr) c 12 '-(tf-p+CT)*^ ^23*23*13 

S 



°12 
C 12 



12 



+ 



P 

Cl3 



( c (<5-p) s 12 s 23 ~~ c p c 12 c 23 s 13) + Cl3 c p ( s (<5-p) s 12 s 23 + s p c 12 c 23 s is) 



c 12 
2 „2 



C 23 l C 12 C 13 



/- ( C (5-cr) C 12 S 23 + C cr C 23 S 12 S 13 ) + C23 C cr ( s (<5-cr) c 12 S 23 S (? C 2Z S \2 S Yi ) 
^23 



c 12 s 13 

2 ^ 

C 23 C 13 S 12 I 



s 13 



dt 167T 2 \ Cl2 



1 2 _ 2 2 

C (p-cr) C 12 S 12 l S 23 C 23 S 13 



) ( C (<5+p-cr) c 12 c (<5-p+cr) s 12j C 23 S 23 S 13 



-12 



'12 



"Cl2 C (p-cr) 



2 _22\_/ 2 1 2 \ 

(p-cr) c 12' 5 12 ^23 — C 23 S 13 J \ S (S+p-a) C 12 ' S (6-p+a) S 12 J C 23 S 23 S 13 



C 



12 



'12 



+ 



{ ( C (5-cr) C 12 S 23 + C <r C 23 S 12 S 13) + C23 C cr ( S (5-cr) C 12 S 23 S o-C 23 S 12 S 13 ^) 
( C (<5-p) S 12 S 23 — C p C 12 C 23 S 13) + Cl3 C p ( S (<5-p) S 12 S 23 + S P C 12 C 23 S ls) 



S 

C23 
Cl3 



c 23 ( C 13 S 12 s 13) 



Si 9S 



12' , 13 



Ci 0C0Q C 



1 2°23°13 
s 13 



(7) 



It is worth mentioning that the analytical results in Eqs. (6) and (7) can also be achieved 
from Eqs. (14)— (19) of Ref. [15] by setting C l K — 1 and y v — over there. One can see 
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that the running effects of three flavor mixing angles and three CP-violating phases are all 
proportional to the factor ?/ 2 /(167r 2 ) = m 2 (l + tan 2 (3)/(1Qtt 2 v 2 ) m 6.6 x 1CT 7 (1 + tan 2 /?) 
in the MSSM. Hence an appreciable value of tan/5 is required, in order to get appreciable 
radiative corrections to relevant physical parameters. 

Next let us consider three rephasing- invariant quantities of V. The first one is the 
Jarlskog parameter J [18], defined through 



7,fc 



(8) 



where the Greek and Latin subscripts run over (e,/x,r) and (1,2,3), respectively. Taking 
account of the parametrization of V in Eq. (1), we have J = Si2 c i2 s 23 c 23 s i3 c i3 s <5- The 
off-diagonal asymmetries of V [19], 



A L = \V e2 \ 2 


- |v^i| 2 = |^ 3 | 2 


-|K 2 | 2 = 


\v Tl 


2 IT/ I 2 
— \V e3 1 , 


A R =\V e2 \ 2 


-|^ 3 | 2 = l^i| 2 


-|K 2 | 2 = 




2 IT/ I 2 


are also rephasing-invariant. 


We obtain A L = 


= s 12 ( c 13 — 


c 23) 


~ C \2 S 23 S \d, 



and A R 
and (7). The result for J is 



-13 \ a 12 



(9) 

2 S 12 C 12 S 23 C 23 S 13 C <5 

— s 2 3 ). The RGEs of J , A L and A R can then be derived from Eqs. (6) 



dj 
dt 



16tt 2 



C 23 C 13 S 23 S 13 



,12 



C S S (p~a) 



-*) S 6 (< 



12 



'12 



C (5+p-o-) C 12 C 23 S 23 S 13 



_ 2 / 2 _ 2 2 \ 

~~ C (<5-p+cr) C 23 S 12 S 23 S 13 C (p-o-) C 12 S 12 \ S 23 C 23 S 13j 



+Cl2 c 23 C 13 S 23 S 13 

+ S (5-p+cr) C 23 S 12 S 23 S 13 

,2 

( ■)■>( 



C S C (p-a) ~ S 6 S (p 



C 23 C 13 S 12 



,13 



p) S 12 S 23 



a) (c? 2 
(p-cr) c 12' 5 'l2 

) 



• s \ ,,_ /' i >- s io (s 



42) 

2 2 
C 23 S 13 



s (5+p-cr) C 12 C 23 S 23 S 13 



C p C 12 C 23 S 13 



" S p C 12 S 12' S 13 l C 23 



( C 23 S 23) 



+ s (5-p) c 23 S 12 s 23 S 13 + 

+Cl3 C 23 C 13 S 12 ( S (<5-p) S 12 

_ 2 2 1 

C (<5-p) c 23 S 12 s 23 S 13 ' 



C 5 S p 



C p S S \ C 13 



'13 



C 12 C 23 S 23 



S 23 + S p C 12 C 23 S 13) 



C C-i nS -1 9S1 



p^l2°12°13 \ C 23 
C 12 C 23 S 23 



c 5 C p + s 5 S p 



S 23) 



12^23^13 



5 23 



S o- c 12 s 12' s 13 \ C 23 
+C23 C 12 C 23 C 13 ( S (<5-cr 



( c 13 ~~ S is) 

( C (<5-cr) C 12 S 23 + C cr C 23 S 12 S 13) { C 8 S a ~ C a S 8 C 13) C 
( C 23 — S 23) + ( S (<5-cr) C 12 + C a S 8 S 12) 



23 S 12 S 23 



(T' 5 5' 5 12 j C 23 S 23 S 13 



-ct) C 12 S 23 



S <t C 23 S 12 S 13 



1 2 \ 2 

C 8 C a ' S 8 S a C 13 J C 23 S 12 S 23 



+ C <r C 12 S 12 S 13 ( C 23 S 23) + ( C (<5-cr) C 12 + S <5 S cr S 12) C 23 S 23 S 13] } 



(10) 



Different from d<5/dt, dj7"/dt does not suffer from any divergence in the #13 — > limit. 
This feature proves that J itself is a well-defined rephasing-invariant quantity, while 5 is 
parametrization-dependent and cannot be well defined when #13 vanishes. The RGEs of A h 
and A R read as follows: 
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df 



2 2 ( 2 

c (<5+p-o-) c 12 c 23' S 23 S 13 ~~ c (5-p+o-) C 23 S 12 s 23 S 13 — c (p-cr) c 12 s 12 \^ s 23 ~~ 



2 2 x 
C 23 S 13, 



167T 2 [Cl2 

_ _l_ ( 2 _ 2 ^ _ f 2 _ 2 2 ^ 

— S <5 S (p-o-) C 23 S 23 S 13 C 5 C (p-<r) C 23 S 23 S 13 \^ C 12 S 12 ^ — C (p-<r) C 12 S 12 ^ S 23 C 23 S 13 J 

i n> I" 2 , 2 / 2 2 2 \ 

~r^S12 [ S (5+p-cr) C 12 C 23 S 23 S 13 + S (5-p+<r) C 23' S 12' S 23 S 13 — S (p-cr) C 12 S 12 1^ S 23 ~~ C 23 S 13 J 



2 2 \ / 2 2 2 

C (p-<r) S 5 C 23 S 23 S 13 + c <5 s (p-<r) c 23 ^ C 12 ~~ s 12 J S 23 S 13 ~~ s (p-cr) c 12 s 12 \^ s 23 ~~ C 23 S 13 

2 C 23 C 13 S 23 



3 13 



p) S 12 S 23 C p C 12 C 23 S 13 



) (<*- 



p) C 23 S 12 + C p C 12 S 23 S 13 



) 



"2Cl3 C 23 C 13 S 2 



? 13 L 23°13' , 23 V S (<5-p) S 12 S 23 + S p C 12 C 23 S 13 ) I S (5-p) C 23 S 12 S p C 12 S 23 S 13 



2 c <t C 23 C 13' S 12 s 13 



C: 



23 



( C (5-cr) C 12 S 23 + C <t C 23 S 12 S 13) 



+ 2C23 S cr C 23 C 13 S 12 S 13 ( S (5-cr) C 12 S 23 S <r C 23 S 12 S 13) } ' 

_ y 2 T cj 3 / 2 C(p- CT )C 12 s 12 



dt 16tt 2 



12 



2 2 
C (<5+p-o-) c 12 c 23 S 23' S 13 + c (5-p+o-) C 23 S 12 s 23 S 13 



2 _ 2 2 

(p-cr) c 12' 5 12 1^23 C 23 S 13, 



12 S (p-cr) C 12' S 12 



• S (<5+p-o-) C 12 C 23 S 23' S 13 



_i 2 / 2 _ 2 2 

+ S (5-P+ct) C 23 S 12 S 23 S 13 S (p-<r) C 12 S 12 l S 23 C 23 S 13 



23° 23 



C: 



13 



( c (5-p) s 12 s 23 + C p C 12 c 23 S 13) ( c (<5-p) c 23 S 12 + C p C 12 s 23 S 13) 



+ 2C13C23S23 ( S (5-p) S 12 S 23 + S p C 12 C 23 S 13 



~ S (5-p) C 23 S 12 + S p C 12 S 23 S 13 



2 C 23 



C23 



( C (<5-cr) c 12 s 23 + C ct C 23 S 12 s 13) 2^23^3 (s(S-a) C 12 s 23 s <r C 23 S 12 s 13) | ■ (H) 



Eqs. (10) and (11) clearly show that the RGE evolution of J7", A L and A R depends on 
the Majorana phases p and a. Hence the radiative corrections to these rephasing-invariant 
parameters would in general be quite different, if neutrinos were Dirac particles instead of 
Majorana particles 3 . 

III. NUMERICAL EXAMPLES AND DISCUSSION 



We proceed to illustrate the radiative generation of three CP-violating phases by taking 
a few typical numerical examples. The eigenvalues of Y\ at Ass are chosen in such a way 
that they can correctly run to their low-energy values [21] 4 . We assume the masses of three 



3 See, e.g., Ref. [20] for a recent analysis of the RGE evolution of Dirac neutrino masses and lepton 
flavor mixing parameters. 

4 A similar treatment is required for the gauge couplings, the quark Yukawa coupling eigenvalues 
and the quark flavor mixing parameters in the full set of RGEs [11-15,22]. 
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light neutrinos to be nearly degenerate and m 1 ~ 0.2 eV, so as to make the RGE running 
effects of relevant physical quantities significant enough. The initial values of k 1 , k 2 and k 3 
can be adjusted via 



m 1 
v 2 sin 2 /3 





+ Am 2 .! 


V 2 


sin 2 (3 




+ Am§! 


V 2 


sin 2 /? 



(12) 



together with a typical input tan/5 = 10, such that the resultant neutrino mass-squared 
differences Am^ = m 2 — m\ and Am^ = m\ — m\ at Aew are consistent with the solar 
and atmospheric neutrino oscillation data [1-5]. We follow a similar strategy to choose the 
initial values of three mixing angles 6*12, #23 and #13, in order to reproduce their low-energy 
values determined or constrained from a global analysis of current experimental data [8]. 
In view of the upper bound #13 < 10°, we shall typically take 6>i 3 = 1°, 3° and 5° in our 
numerical calculations. We allow one of three CP-violating phases (5, p, a) to vanish at Ass 
and examine whether it can run to 90° at Aew by choosing the initial values of the other 
two phase parameters properly. 



A. Radiative Generation of S = 90° 

First of all, let us look at the one-loop RGE evolution of 5 in the m 1 < m 2 < m 3 case 
(i.e., Arajj > 0). The input and output values of relevant physical parameters are listed in 
Table I. One can see that 5 = 90° at Aew can be radiatively generated from 5 = 0° at Ass, 
if 6*13 = 1°, p = 4.0° and a = —57.5° are input. Changing the initial value of 6*13 to 3° or 5° 
but fixing the input values of the other quantities, we find that only 5 = 41.8° or 5 = 35.8° 
can be obtained at Aew- While the results of m 1 , Am 2 1; #i 2 and 6*23 at Aew are rather 
stable against the change of #13 from 1° to 5° at A S s, the result of Am^ becomes smaller 
and less favored. 

The running behaviors of 5, p and a are explicitly shown in FIG. 1, in which those of J~, 
A L and ^4 R are also illustrated. Some comments are in order. 

• The one-loop RGE running behaviors of three CP- violating phases are quite similar in 
the chosen parameter space. It is easy to understand this feature from Eq. (7): given 
the conditions that holds and 6> 13 is small, Eq. (7) can be simplified to 



dS 
dt 
dp 
dt 
da 
dt 



16tt 2 

Jl_ 

16tt 2 

16tt 2 



c (p-a) S (p-a) 2 



C: 



'23 



+ 



C (S-P) S P 



C (8-a) S a 



12 



C: 



13 



23 



C 12 S 12 C 23 S 23 



'13 



C (p-a) S (p-a) 2 2 , 
* 12-^23 



C (S-p) S p C (S-a) S a\ C l2 S 12 C 23 S 
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C 12 6 23 ~r 



C: 
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C: 



13 



C: 



23 



'13 
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Thus the trend of p or o in the RGE evolution is very similar to that of 5. 

• The Jarlskog parameter is sensitive to both 5 and # 13 . Given the initial condition 
#13 = 1°, 5 = 90° at A E w can be achieved from 5 = 0° at A ss . It turns out that 
J 0.26% can be radiatively generated in this specific case. When 9\ 3 = 5° is input 
at Ass, one may arrive at J ~ 1% at Aewj a value which might be detectable in the 
future long-baseline neutrino oscillation experiments. 

• The off-diagonal asymmetries A L and ^4 R are both non-vanishing in our numerical 
example. While A L is sensitive to the input of #13, A R is not. The reason for the latter 
feature is simply that ^4 R = c\ 3 (s\ 2 — s\ 3 ) holds. 

It is worth mentioning that the unitarity of V requires that its six unitarity triangles in the 
complex plane have the same area, amounting to J /2 [22]. Thus the radiative generation 
of 5 or J geometrically implies the radiative generation of every leptonic unitarity triangle; 
i.e., three overlapped lines (sides) at Ass can evolve into a triangle at Aew- 

Now let us examine whether the radiative generation of 5 = 90° at Aew can be achieved 
from other initial values of p and a, when 9 13 = 1° holds and m 1 , Am^, Am^, 12 and 9 23 
take the same input values as before at Ass ( see Table I). We find out two new numerical 
examples with (p, a) = (0°, —59.9°) and (10°, —57°), respectively, as shown in FIG. 2. Note 
that both S and p (or a) can be radiatively generated from a 7^ 0° (or p 7^ 0°) at the seesaw 
scale 5 . Our results hint at the existence of strong parameter degeneracy in obtaining 5 = 90° 
at Aew from 5 = 0° at Ass- To resolve this problem is certainly a big challenge in model 
building, unless two Majorana CP- violating phases could separately be measured at low 
energies. 

Finally we demonstrate that it is also possible to radiatively generate 5 = 90° at A E w 
from 5 = 0° at Ass in the Am^ < case. Our numerical results are shown in FIG. 3, where 
the initial values of relevant parameters are listed in Table II. One can see that the one-loop 
RGE running behaviors of three CP- violating phases in FIG. 3 are very similar to those in 
FIG. 1, although the initial conditions in these two cases are quite different. The reason for 
this similarity has actually been reflected by Eq. (13), where the sign of Am^ does not play 
a role in the leading-order approximation. 

B. Radiative Generation of p = 90° 

We continue to take a look at the one-loop RGE evolution of p in the m 1 < m 2 < m 3 
case (i.e., Am^ > 0). The input and output values of relevant physical parameters are 
listed in Table III. One can see that there is no difficulty to radiatively generate p = 90° at 
A E w from p = 0° at A S s, provided 9 13 = 1°, 5 = 0° and a = —67.7° are input. Allowing the 



5 The possibility to simultaneously generate p and a from 8 7^ 0° at Ass V18k the RGE evolution 
is in general expected to be strongly suppressed, because the leading terms of dp/dt and da/dt 
in Eq. (7) vanish for p = a = 0° at Ass- Only in the 9\ 3 — > limit, the running effects of three 
CP-violating phases could become significant [15]. 



9 



initial value of #13 to change to 3° or 5° but fixing the input values of the other quantities, 
we find that only p = 17.6° or p = 12.1° can be obtained at A E w- Again the results of m 1 , 
Ara^, 612 and #23 at A EW are quite stable against the change of #13 from 1° to 5° at Agg, 
but the result of Am^ becomes smaller. 

The RGE running behaviors of three CP- violating phases (5, p, a) and three rephasing- 
invariant parameters (J7, A L , A R ) are explicitly shown in FIG. 4. Comparing this figure 
with FIG. 1, one may observe much similarity between them. There are two reasons for this 
similarity: first, the initial conditions in these two cases are not very different; second, the 
evolution of 5, p and a is dominated by Eq. (13) in both cases. 

C. Radiative Generation of a = 90° 

Finally let us look at the RGE evolution of a in the m 1 < m 2 < m 3 case (i.e., Am^ > 0). 
The input and output values of relevant physical parameters are listed in Table IV. One can 
see that there is no difficulty to radiatively generate a = 90° at Aew from a = 0° at Ass, if 
#13 = 1°, 5 = 119.7° and p = 60.8° are input. Allowing the initial value of #13 to change to 
3° or 5° but fixing the input values of the other quantities, we find that only a = 34.4° or 
a = 29.8° can be obtained at A EW . The results of m 1 , Aro^, 9 12 and #23 at A EW are very 
stable against the change of 6*13 from 1° to 5° at A S s, but the result of Ara^ becomes larger 
and less favored. 

The RGE running behaviors of 5, p and a are shown in FIG. 5, in which those of J7, A L 
and A K are also illustrated. Again the evolution of three CP-violating phases is dominated 
by Eq. (13). Because 5 evolves from the second quadrant to the third one during the RGE 
running, a flip of the sign of J appears around S = 180°. The signs of A L and A R keep 
unchanged from Ass to A E w, implying that the geometrical structure of V does not change 
in a significant way. 

In all the cases discussed above, (m) ee ~ m 1 holds as a consequence of the approximate 
mass degeneracy of three light neutrinos. Indeed, the smallness of 9 13 allows us to obtain 
an approximate expression of (m) ee : 

(m) ee ps m 1 yl — sin 2 2# 12 sin 2 (p — a) G [m 1 cos 26> 12 , wij] . (14) 

Given m l m 0.2 eV and 9i 2 ~ 33° at low-energy scales, {m) ee turns out to lie in the range 
0.08 eV < (m) ee < 0.2 eV, a result consistent with the present experimental upper limit 
(m) ee < 0.38 eV [8]. 

IV. SUMMARY 

Taking account of a very useful parametrization of the 3x3 MNS matrix V, we have 
derived the one-loop RGEs for its three mixing angles (#12, #23, #13), three CP-violating 
phases (5, p, a) and three rephasing-invariant quantities (J7, A L , A-^). Particular attention 
has been paid to the radiative generation of leptonic CP violation, because the Dirac phase S 
and the Majorana phases p and a are entangled with one another in the RGE evolution from 
the seesaw scale A S s to the electroweak scale A EW . We have shown that 5 = 90° at A EW can 
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be radiatively generated from 8 = 0° at Ass i n the minimal supersymmetric standard model, 
provided the light neutrino masses are nearly degenerate and the mixing angle 6 13 is of 0(1°) 
or smaller. As for p and a, it is also possible to radiatively generate p = 90° or a = 90° at 
A EW from p = 0° or a = 0° at A S s- An interesting feature of our numerical analysis is that 
the one-loop RGE running behaviors of three CP-violating phases are quite similar in the 
chosen parameter space, no matter whether the sign of the neutrino mass-squared difference 
Am| x is positive or negative. As an important by-product, the geometrical structure of 
V (i.e., its off-diagonal asymmetries and the area of its unitarity triangles in the complex 
plane) against radiative corrections has also been discussed. 

At this point, it is worthwhile to summarize some remarkable differences between our 
analytical and numerical results and those obtained in Refs. [12,13]: 

(1) The phase convention of V used in Refs. [12,13] leads to the apparent dependence of 
(ra) ee (the effective mass of the neutrino less double-beta decay) on 5, thus it is ill to refer 
to 5 as the "Dirac" CP-violating phase. Such an ambiguity has been avoided in our phase 
convention of V, which forbids 5 to enter the expression of (m) ee . Our Majorana phases (p 
and a) are related to the ones defined by Casas et al [12] (0 and <fi') or by Antusch et al [13] 
((p l and ip 2 ) as follows: tp 1 — (f> — 2 (5 — p) and tp 2 — (f)' — 2 (5 — a). On the other hand, 
our Dirac phase 5 is indistinguishable from their phase parameter 5 in the description of 
leptonic CP violation in neutrino oscillations. 

(2) Our RGEs for three mixing angles and three CP-violating phases are apparently 
different from those given in Refs. [12,13], just because we have taken a different and more 
instructive phase convention for V. In particular, the leading terms of dp/dt and da/dt (i.e., 
the terms proportional to I/C12 an d V s i3) are very similar to those of d8/dt, as one can see 
from Eq. (7) or Eq. (13). This result implies that three CP- violating phases have similar 
RGE running behaviors in our phase convention. In contrast, the RGE running effect on tp l 
(or 0) and (p 2 (or <$') is much weaker than that on 5 in Refs. [12,13]. 

(3) The one-loop RGEs for J~, Ah and A-r given in Eqs. (10) and (11) are new results, 
although they can be derived from Eqs. (6) and (7) in a straightforward way. These three 
quantities are actually rephasing-invariant or parametrization-independent. Hence we can 
get some more generic feeling about radiative corrections to the leptonic CP-violating effect 
and the geometrical structure of V from the RGE evolution of J, Ah and ^4 R . This kind of 
study was not done in Refs. [12,13] or elsewhere. 

(4) Different from those previous works, we have concentrated on the novel possibilities 
to radiatively generate 5 = 90°, p = 90° or a = 90° at A EW from 5 = 0°, p = 0° or 
a = 0° at A ss in our numerical exercises. Because of J oc sin 5, 5 ~ 90° at low energies 
is a necessary condition to achieve sufficiently large CP- and T-violating effects in neutrino 
oscillations. On the other hand, p ~ 90° or a ~ 90° at low energies may result in a kind 
of large cancellation in (m) ee , implying the possible suppression of this unique experimental 
observable to identify the Majorana nature of massive neutrinos. 

It is also worth remarking that our analysis is essentially independent of the specific 
textures of lepton Yukawa coupling matrices, thus it can be applied to the concrete work 
of model building. Since the elegant thermal leptogenesis mechanism is usually expected 
to work at the seesaw scale, a study of its consequences at low-energy scales is available 
by means of the one-loop RGEs that we have obtained. In other words, the RGEs may 
serve as a useful bridge to establish a kind of connection between the phenomena of CP 
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violation at low- and high-energy scales. It would be extremely interesting, in our opinion, 
if the phase parameter governing the strength of CP violation in a long-baseline neutrino 
oscillation experiment could be radiatively generated from those CP-violating phases which 
control the matter-antimatter asymmetry of our universe at the seesaw scale. 

One of us (Z.Z.X.) is grateful to S. Zhou for very helpful discussions. This work was 
supported in part by the National Nature Science Foundation of China. 
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TABLES 



TABLE I. Radiative generation of 5 = 90° at the electroweak scale Aew from 5 = 0° at the 
seesaw scale Ass in the MSSM with tan (5 = 10 and Am^ > 0. 



Parameter 


Input (Ag S ) 




Output (Aew) 








0n = 1° 


9i* = 3° 

w lo ,f 


#iq = 5° 

w lo J 


m 1 (eV) 


0.241 


0.20 


0.20 


0.20 


Am^lO- 5 eV 2 ) 


20.4 


7.79 


7.17 


6.56 


Aml^lO- 3 eV 2 ) 


3.32 


2.20 


2.20 


2.20 


#12 


24.1° 


33.0° 


33.0° 


33.1° 




43.9° 


45.1° 


45.0° 


45.0° 


#13 


l /375 


0.65° 


2.46° 


4.52° 


5 


0° 


90.0° 


41.8° 


35.8° 


P 


4.0° 


72.2° 


23.8° 


17.6° 


a 


-57.5° 


26.3° 


-22.0° 


-28.1° 


TABLE II. Radiative generation of 5 = 


90° at the electroweak scale Aew from 6 


= 0° at the 


seesaw scale Ass 


in the MSSM with tan/3 = 


10 and Am!^ < 0. 






Parameter 


Input (A ss ) 




Output (Aew) 








#13 = 1° 


#13 = 3° 


#13 = 5° 


m 1 (eV) 


0.241 


0.20 


0.20 


0.20 


Am^RT 5 eV 2 ) 


20.8 


7.96 


7.37 


6.79 


Amf^Kr 3 eV 2 ) 


-3.08 


-2.20 


-2.21 


-2.21 


#12 


25.1° 


33.3° 


33.5° 


33.7° 


#23 


46.7° 


45.2° 


45.2° 


45.2° 


#13 


l73 /5 


0.68° 


2.38° 


4.34° 


5 


0° 


90.0° 


43.1° 


36.8° 


P 


-81.9° 


-13.6° 


-60.7° 


-67.2° 


a 


34.3° 


117.7° 


70.6° 


64.1° 
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TABLE III. Radiative generation of p = 90° at the electroweak scale Aew from p = 0° at the 
seesaw scale Ass in the MSSM with tan (3 = 10 and Am^ > 0. 



Parameter 


Input (A ss ) 




Output (Aew) 








#13 = 1° 


#13 = 3° 


#13 = 5° 


m 1 (eV) 


0.241 


0.20 


0.20 


0.20 


Am^ltr 5 eV 2 ) 


20.4 


8.54 


7.90 


7.27 


Ami^KT 3 eV 2 ) 


3.32 


2.21 


2.20 


2.20 


#12 


27.6° 


33.1° 


33.2° 


33.3° 


#23 


43.9° 


44.8° 


44.8° 


44.8° 


#13 


l73 /5 


0.43° 


2.17° 


4.24° 


5 


0° 


107.6° 


35.4° 


30.2° 


P 


0° 


90.2° 


17.6° 


10 10 

12.1 


a 


-67.7° 


34.1° 


-38.3° 


-44.6° 


TABLE IV. Radiative generation of a = 


= 90° at the electroweak scale Aew from a 


— 0° at the 


seesaw scale Ass 


in the MSSM with tan (3 = 


-- 10 and Am^ > 0. 






Parameter 


Input (A ss ) 




Output (Aew) 








#13 = 1° 


#13 = 3° 


#13 = 5° 


m 1 (eV) 


0.241 


0.20 


0.20 


0.20 


Am^RT 5 eV 2 ) 


20.3 


8.24 


8.85 


9.50 


Am^KT 3 eV 2 ) 


3.32 


2.21 


2.21 


2.21 


#12 


24.1° 


33.2° 


34.2° 


35.2° 


#23 


43.9° 


44.6° 


44.6° 


44.6° 


#13 


r/3 /5 


0.51° 


2.27° 


4.29° 


5 


119.7° 


216.0° 


161.0° 


157.1° 


P 


60.8° 


135.3° 


78.8° 


73.6° 


a 


0° 


90.0° 


34.4° 


29.8° 
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FIGURES 
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FIG. 1. The RGE running behaviors of three CP-violating phases and three rephasing-invariant 
quantities of V from Agg to Aew m the MSSM, where the input values of relevant parameters can 
be found from Table I. 
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FIG. 2. The radiative generation of 5 = 90° at Aew from #13 = 1° and different values of 
and a at Ass- The input values of other parameters can be found from Table I. 
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FIG. 3. The RGE running behaviors of three CP-violating phases and three rephasing-invariant 
quantities of V from Ags to Aew i n the MSSM, where the input values of relevant parameters can 
be found from Table II. 
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FIG. 4. The RGE running behaviors of three CP-violating phases and three rephasing-invariant 
quantities of V from Ags to Aew i n the MSSM, where the input values of relevant parameters can 
be found from Table III. 
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FIG. 5. The RGE running behaviors of three CP-violating phases and three rephasing-invariant 
quantities of V from Ags to Aew i n the MSSM, where the input values of relevant parameters can 
be found from Table IV. 
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